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ABSTRACT 

Detailed calculations are described of the rate of heating of 
the ambient electrons arising from solar ultraviolet radiation and 
it is shown that the resulting difference between the electron tem- 
perature and the gas temperature at noon in a quiet ionosphere 
may reach a maximum value lying between 800° K and 1500° K at 
an altitude near 2 00 km but that the difference vanishes below 
about 120 km and above about 400 km. This is in accord with the 
analysis of the observational data by Bauer and Bourdeau which 
can therefore be explained as due to the direct action of solar 
ultraviolet radiation and it is unnecessary to postulate the exist- 
ence of electric fields. 
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ELECTRON TEMPERATURES IN THE IONOSPHERE 


A. Dalgarno, 

M. B. McElroy* and R. J. Moffett* 

1. Introduction 

Hanson and Johnson (1961) have suggested that energetic photo- 
electrons produced by the incident solar ultraviolet radiation may give 
rise to an electron temperature Te in the daytime ionosphere which is 
higher than the neutral particle and positive ion temperature Tg and 
experimental evidence supporting this view has been obtained recently 
by Spencer, Brace and Carignan (1962) who find that Te is much greater 
than Tg at altitudes between about 100 and 400 km, the maximum value 
of Te lying between 2400°K and 3000°K. However, because the solar 
photon flux adopted by Hanson and Johnson is much less than that now 
observed (cf. Watanabe and Hinteregger 1962), because an important 
process--the electron impact excitation of vibrational levels of molec- 
ular nitrogen — was omitted from consideration, and because the significance 
of metastable energies as a source of heating was not investigated, the 
conclusions of Hanson and Johnson have only qualitative significance and 
a re-analysis is desirable. 


*Department of Applied Mathematics, The Queen's University of Belfast, 
Belfast, N. Ireland. 
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2. Production of Photoelec trons 


2.1 Model Atmospheres 

As model atmospheres we have selected two from those derived by 
Bates (1959) and Bates and Patterson (1961), one of which corresponds 
to a temperature of 1000°K at the base of the exosphere and the other 
to a temperature of 2000°K, The two models probably span the possible 
atmospheres. The associated number densities and temperatures are 
reproduced in Table 1. 

2.2 Solar Flux and Photoionization Rates 

For the incident flux of solar photons we have adopted the values 
quoted by Watanabe and Hinteregger (1962) and they are reproduced in 
Table 2. Using a model atmosphere which is roughly a mean of the two 
given in .Table 1, Watanabe and Hinteregger have computed the rate of 
production of ionization corresponding to this solar flux. We have 
repeated their calculations for the model atmospheres in Table 1 and a 
comparison of the noon ionization rates is given in Figure 1. The rela- 

- j— ■'P -j- 

tive production of 0 , 0^ and N ? ions are of interest and these are 
shown in Figures 2 and 3. The accuracy of the absorption and photo- 
ionization cross sections used is very uncertain and Figures 1 and 2 
should be regarded as merely representative of the actual ionization 
rates . 
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Atmospheric Temperatures and Number Densities 




y“N 

y-N 

y — \ 

y — s 

y-^ 

y~% 

o 

rH 





00 

CO 

00 

CO 

CO 

H 

iH 


/*N 

v-y 

v-y 

V-/ 

v_y 

^y 

\ y 

N_y 

s»y 

vy 


CM 

CM 

CO 

o 

o 

O 

rH 

CO 

vO 

o 



iH 

uo 

st 


CO 

CO 

CM 

ON 

st 


c 

CO 

v£> 

H 

oo 

rs 

rH 

vC 

OO 

St 











y-N 



y~s 

/"“V 


y-N 

y—s 

y—s 

y—s 

y-N 

O 



vO 

vC 


rs 


oo 

CO 

ON 

I J 



N_y 

v_y 

Vy/ 

v_y 

v^ 

'w y 

H' 

v-y 

v — y 


o 

2 

CM 

CO 

CM 


rH 

rs 

in 

00 

o 


o o 

O 

O 

vD 

r— 1 

s- 

CO 

rH 

CM 

00 

00 


© c 

CM 

CM 

st 

rH 

CM 

is. 

CM 

00 

vO 

co 


ti 









y-N 



y-N 


y— N 

y— v 

y—s 

y — n 

y-N 

y-N 

o 



rs 

00 

00 

00 

00 

00 

co 

ON 

rH 


8 o 

Nw^ 

v-y 

v^ 

V_y 

v_y 

^ — y 

v-y 

N_y 

N-/ 


w v_y 

<r 

ON 

o 

o 

ON 

rs 

CM 

co 

o 


H C 

rH 

o 

r-» 

rs 

St 

O 

rs 

ON 

st 




iH 

rH 

CM 

st 

co 

rH 

U0 

oo 


y-N 












O 

o 

o 

n 

O 

o 

o 

o 

o 


o 

o 

is. 

st 

on 

iH 

o 

on 

00 

00 



o 

ON 

C A 

CO 

00 

vC 

oo 

oo 

oo 


to 

#v 

o 

r\ 

«v 

•s 

rv 

*> 




H j 

i 

CM 

iH 

rH 

rH 

rH 

rH 

rH 





/-N 

y—s 

y— v 

y-N 

y-N 

y“N 

y-N 

y-N 

y-N 



vO 

vC 

vC 

vO 

o 

vC 

U0 

uo 

in 



N_y 

V-/ 

V 

v_ y 

v^y 

s-y 

v-y 

v-y 

' — y 


O 

m 

St 

st 

00 

vC 


O 

O 

o 


c 

o 

<M 

St 

vC 

vD 

rH 

CM 

vD 

CM 



rH 

rH 

iH 

rH 

iH 

rH 

rs 

St 

00 










y-N 

y-N 



y-\ 

/-N 

y — v 

/— n 

y—s 

y-s 

y— N 

o 

rH 



UO 

vC 


p-n 

00 

ON 

ON 

iH 

rH 


/-N 

S-' 

v_y 

V^ 

V— y 

vy 

s— y 

v — ' 

v— y 

v-y 


CM 

VO 

co 

vC 

St 

00 

C7\ 

ON 

O 

O 


523 

00 

VC 

rH 

rH 

00 

O 

uo 

CM 

st 


N—' 

• 

• 

• 

. 

. 

. 


. 



c 

vO 

CM 

rH 

U0 

CM 

iH 

uo 

st 

St 











y-N 



y-\ 

/-N 

y-N 

y — n 

y--\ 

y“N 

y— N 

y-N 

o 


<-s 

st 

U0 

U0 

vC 


oo 

CO 

ON 

rH 


! 'w' 

' — y 

v^ 

V 

v_y 

v 

N_y 

' — y 

v-y 


O CM 

st 

«H 

CM 

St 

CM 

vC 

OO 

00 

o 


o o 

O W 

rH 

rH 

CO 

CM 

00 

O 

00 

vD 

CO 


o a 

r> 

rH 

CM 

iH 

U0 

O0 

iH 

>H 

vO 

vD 

00 

OO 

1 1 

II 

/H 

y*N 

y-N 

y ~ n 

y—\ 

y-N 

y-N 

y-N 

y-N 

o 

B 

M 



O'* 


CO 

co 

00 

ON 

ON 

rH 

w 

/£N y-N 

V-y 

v«/ 

V/ 

N_y 

v^/ 

v-y 

v-y 

v-y 

vy 


8 o 

oo 

VC 

st 


CM 

CM 

ON 

ON 

O 

i 

o 

W V 

H C 

O 

CO 

st 

CM 

O 

00 

00 

UO 

st 

rH 

iH 

CM 

uo 

tH 

00 

rs 

rH 

vO 

CO 

x 











oo 











o 

/-N 

O 

O 

o 

o 

CO 

O 

CM 

oo 

o 

rH 

o 

o 

o 

o 

o 

on 

co 

uo 

vO 

00 

II 

V-y 

o 

O 

o 

o 

co 

ON 

ON 

is 

00 

to 

H 

rH 


rH 

rH 






y-N 

rs 











vy 











OO 

a 










O 

TD 

4J £ 

o 

o 

o 

o 

o 

o 

o 

O 

o 

rH 

•H }2 

o 

UO 

o 

U0 

o 

in 

o 

uo 

CM 


4J W 
pH 
< 

UO 

st 

sr 

00 

oo 

CM 

CM 

i — l 

i — i 

+ 


3 



Table 2 


Solar Spectral Intensities at Normal Incidence 


X(in A) 

Flux (in lO^cm ^sec 

X(m £) 

1A 8 -2 -1. 
Flux (xn 10 cm sec ) 

1025.7 

26 

500 - 600 

30 

1000 - 1027 

15 

400 - 500 

24 

989.8 

5 

303.8 

43 

977.0 

30 

300 - 400 

29 

972.5 

10 

230 - 300 

31 

94-9.7 

5 

170 - 230 

33 

911 - 1000 

37 

110 - 170 

3.5 

850 - 911 

95 

80 - 110 

2.4 

796 - 850 

25 

60 - 80 

1.8 

700 - 796 

50 
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1.5 

600 - 700 

47 
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0.12 

584.3 

29 

10 - 20 

0.02 
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Figure 2. Noon photoionization rates for T (co) = 1000 
for 0 + , 0 + and N + . 





Figure 3. Noon photoionization rates for 
for 0 + . 0 + and K’ + . 





2.3 Equilibrium Electron Densities 


Since the rates of the processes responsible for the disappearance 
of free electrons are still very uncertain, we have not attempted to 
calculate the equilibrium distributions of the electron density as 
functions of altitude for the assumed models but we have instead adopted 
the values given by Ratcliffe (1960) and listed in Table 1. 

2.4 Energy Distribution of Photoelectrons 

The prediction of the energy distribution of the photoelectrons is 
complicated by the fact that the available cross section data refer to 
a multiplicity of ionizing processes whereas the energy distribution 
depends upon the individual processes. Thus, the spectral heads for 
the ejection of an outer shell electron from atomic oxygen for the 
transitions 

0(ls 2 2s 2 2p 4 ) 3 P + hv - 0 + (ls 2 2s 2 2p 3 ) 4 S + e (1) 

0(ls 2 2s 2 2p 4 ) 3 P + hv - 0 + (ls 2 2s 2 2p 3 ) 2 D + e (2) 

0(ls 2 2s 2 2p 4 ) 3 P + hv - 0 + (ls 2 2s 2 2p 3 ) 2 P + e (3) 

are located at 910^, 732^ and 663i(, respectively, so that a photon of 

o 

wavelength, say 500A, absorbed by atomic oxygen will produce an electron 
of energy 11.2eV, 7.9eV or 6.1eV depending upon which transition is 
involved. The internal energies of the me tas table ions may still be 
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I 


2 

available as heating. Thus, the radiative lifetime of the state 

3 2 4 

is 5 x 10 sec and of the state is 2 x 10 sec; and the rate 

coefficient for the superelastic collision 

e + 0 + ( 2 D) -* e + 0 + ( 4 S) (4) 

-8 3 -1 

is about 3 x 10 cm sec (Seaton and Osterbrock 1957). The meta- 

2 

stable energy of the D state may therefore be transferred directly 

to the ambient electron gas. Alternatively, the metastable ion may 

undergo ion-atom interchange with the atmospheric molecules, its 

energy thereby appearing as thermal energy of the heavy particle gas. 

2 2 

The radiative lifetime of the P^/2 state ^ sec anc ^ t * ie ^ 3/2 
state is 4 sec and the rate coefficients for the superelastic colli- 
sions 


e + 0 x ( P) 


e + 0 + ( 2 ?) 


e + 0 + ( 2 D) 


e + 0 + ( 4 S) 


(5) 

( 6 ) 


-7 3 -I -8 3 -I 

are, respectively, about 1 x 10 cm sec and 2 x 10 cm sec (Seaton 

and Osterbrock 1957). Except at low altitudes where ion-atom interchange 

2 

may occur more rapidly, the P state is usually deactivated by spontaneous 
emission, the more probable process being 

0 + ( 2 P) - 0 + ( 2 D) + hv (7) 
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yielding a photon of wavelength 7330 S and leaving a metastable 
state, the energy of which may be transferred to the ambient electrons 
or to the heavy particles. 

At 435 the photoionizing transition 

0(ls 2 2s 2 2p 4 ) 3 P + hv-f 0 + (ls 2 2s2p 4 ) 4 P + e (8) 

and at 310 S, the photoionizing transition 

0(ls 2 2s 2 2p 4 ) 3 P + hv -*■ 0 + (ls 2 2s2p 4 ) 2 P + e, (9) 

which involve the ejection of an inner shell electron, can occur 
(Dalgarno and Parkinson 1960). The excited state produced by (8) will 
subsequently decay according to 

4 * 2244 4 - 2234 

0 (Is 2s p ) P -+ 0 (Is 2s 2p ) S + hv (10) 

with the emission of a photon of wavelength 833 £ which can be 
absorbed by atomic oxygen yielding an electron with an energy of 1.2 eV 
or by molecular oxygen yielding an electron of energy 2.7 eV. The 
excited state produced by (9) will decay according to 

4 - 2 42 4 - 2222 

0 (Is 2s2p )T> 0 (Is 2s 2pTD + hv (11) 

or according to 

0 + (ls 2 2s2p 4 ) 2 P 0 + (ls 2 2s 2 2p 3 ) 2 P + hv . (12) 
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The photons emitted in (11) and (12) have wavelengths of 581 £ and 
537 $, .respectively, and can produce further ionization. As before, 
the metastable states in (11) and (12) may be deactivated by super- 
elastic collisions. 

Radiation of wavelengths between 310 & and the K-shell absorption 
edge at about 23 & can produce electrons with energies ranging from 
zero to 540 eV. These high energy electrons will in turn produce 
further electrons by impact ionization. For initial energies greater 
than perhaps 70 eV the mean number of ion pairs produced can be obtained 
by dividing the initial energy by 34 eV but below 70 eV the mean energy 
required to produce an ion pair must rise rapidly as an increasing frac- 
tion of the energy is expended in excitation processes (cf. Da Igarno and 
Griffing 1958). In computing the ionization rates of Figure 1, we have 
arbitrarily assumed that for electron, energies between 20 and 70 eV, one 
ion pair is produced. 

We shall not discuss the absorption of radiation shorter than 23 $ 
in detail. For our purposes it suffices to note that the solar flux of 

photons with wavelengths shorter than 30 & is of the order of 0.1 erg 

-2 -1 

cm sec and that it is largely deposited in the E-region of the 
ionosphere. 

The photoionization pattern for the atmospheric molecules is 
basically similar to but, because of the vibrational structure, more 
complex than that of atomic oxygen. 
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The first ionization potential of molecular nitrogen is equivalent 

o 

in energy to a wavelength of 796 A and corresponds to the removal of a 


0 g2p electron from the ground state configuration 


N(KK(d 2 s) 2 (0 2s) 2 (7T 2p) 2 (d 2p 2 )) 
2\ g u u g / 


g 


+ 2 + 
leaving an molecule in its ground electronic 2 state. At a 

2m g 

slightly shorter wavelength, a TT^p electron can be ejected, yielding 

I 2 p 

N. in an excited TT state, and at 661 A, an inner shell d (2s) can 
2 u u 

+ 2 + 

be ejected, yielding in an excited 2^ state. Presumably both 

2 2 + 

the and 2^ states are deactivated by spontaneous emission, the 
2 + 

2^ giving rise to the first negative band system. There will be 
further discontinuities in the photoionization cross sections 


associated with the removal of a d (2s) electron and yielding N in 

g ^ 

2 + 

an excited 2 state. None of these excited states is metastable 


g 

and there seems to be no possibility of a direct transfer of their 
energy to the ambient electrons. They are of interest, however, 
within the context of the heat balance of the high atmosphere for 
they provide a mechanism by which some of the incident solar energy 
is converted directly into radiation and so rendered unavailable for 
heating. 

The first ionization potential of molecular oxygen is equivalent 
in energy to a wavelength of 1016 % and corresponds to the removal of 
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a tt 2p electron from the 0 ground state configuration 
g 2 

0 2 (^(d g 2s) 2 (O u 2s) 2 (O g 2p) 2 (7r u 2p) 4 (- g 2p) 2 ") 3 2‘ 


+ o 

leaving an 0^ molecule in its ground electronic tt ^ state. At a 

wavelength of perhaps 767 £, a TT^p electron can be ejected, yielding 
+ 4 

an 0- molecule in the excited tt state, and at a wavelength of 
2 u 

729 it can yield an CL + molecule in the excited V state. The 

2 u 

d 2p electron can be ejected at a wavelength of perhaps 681 yieid- 
S 

ing 0 + in the ^2 state and there will be other discontinuities 
2 g 

at shorter wavelengths associated with the removal of d 2s and d 2s 

u g 

4-2 2 2 

electrons, yielding further 0. in tt and tt states. The tt and 

2 u g u 

4 - 

2 states are the upper states of the first and second negative 
S 

band systems, the latter of which terminates in the ^77 state. The 

u 

4 

77^ state is metastable and its energy of about 4eV may be transferred 
directly to the ambient electron gas by superelastic collisions. 

The transformation of this essentially qualitative discussion 
into a quantitative prediction of the energy distribution of the 
photoelectrons must be somewhat arbitrary. For the photoionization 
cross sections, we have used the values listed by Watanabe and 
Hinteregger (1962) and we have assumed that whenever the photon 
energy is such that a multiplicity of transitions can occur the 
probability of a particular transition is proportional to the 
statistical weight of the product positive ion state. The transitions 
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which have been taken into account are given in Table 3. 

The resulting distributions with altitude of the rates of 
production of photoelectrons of various energies arising from the 
transition (1) are illustrated in Figures 4 and 5 for the two model 
atmospheres. Similar curves may be constructed for the other 
photoionizing processes. The total rates of deposition of photo- 
electron kinetic energy are given in Figures 6 and 7 which also show the 

distributions with altitude of the rates of deposition of energy 

+ 2 +4 

retained initially in the metastable states 0 ( D) and 0^ ( 7r u ) • Part 
of the energy will subsequently be converted into radiation and the 
remainder into heat, the division depending upon the collision 
processes by which the fast electrons are slowed down. 
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Table 3 

Transitions Taken Into Account in Deriving the Energy 
Distribution of Photoelectrons 


Atomic Oxygen Transitions 


Electron removed 
2p 
2p 
2p 
2s 
2s 


Resulting state 


4 

2 

2 

4 

2 


S 

D 

P 

P 

P 


of 


+ 

0 


Threshold energy in eV 

13.6 

16.9 

18.7 

28.5 

40.0 


Molecular Oxygen Transitions 


+ 


Electron removed 

Resulting state 

o 

Mh 

0 

1 ■ 

Threshold energy 

in eV 

T r 2p 
g 

4 

g 

a 

12. 1 


TT 2p 
u 

u 


16.2 


TF 2p 
U 

? 

~TT 

u 


17 .0 


0 2p 

g 

V 

g 


18.2 


d 2s 
u 

2 ir 

u 


~ 28 


d 2s 
g 

2 rr 

g 


~ 40 


Molecular nitrogen 

transitions 




Electron removed 

Resulting state 

Of N + 

Threshold energy 

in eV 

d 2p 

g 

2 2 + 

g 

2 

15.6 


rr 2p 

u r 

2 t t 

u 


16.7 


d 2s 
u 

V 

u 


18.8 


d 2s 

a 

to 

V 

g 


~ 35 
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Figure 5. Distribution of photoelectron sproduced by process (1) 
for T (») = 2000%. 
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Figure 7. Total rate of deposition of photoelectron kinetic 
energy (curve a) and of metastable, energy (curve 
for T («) « 2000°K. 



3. Energy Loss Processes 


3.1 Excitation and Ionization 

The rate of energy loss per cm of path by optically allowed 

excitation and ionization processes for a high energy electron moving 

-3 . 

through air with a number density n cm is given by the Betne 
formula (ignoring relativistic effects) 

dE 1.87x10 f) /p ^ Z' 1 *5 \ 

— = - n (E/oU) ev cm Uv 

dx E 

where E is the energy in eV (cf. Dalgarno 1961). The Bethe formula 
(13) remains valid to well within an order of magnitude down to 
energies of five hundred electron volts. To obtain the energy loss 
rate at lower energies we assume that for incident energies greater 
than 50 eV, the efficiencies of ionization processes and excitation 
processes are equal and that the mean energy loss for each ionizing 
collision is 20 eV. The rate of energy loss can then be obtained 
using the measured values of the cross sections for the electron 
impact ionization of the atmospheric gases (Tate and Smith 1932, Fite 
and Brackman 1959, Rothe, Marino, Neynaber and Trujillo 1962). At 
500 eV, this procedure yields an energy loss rate in air of 
6 x 10 ^ n eV cm ^ in harmony with that derived from (13). For 
energies between 50 eV and 20 eV we assume that the efficiency of 
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excitation in slowing down the electron is twice that of ionization 
and that the mean energy lost in each ionizing process is 15 eV. 

These procedures are somewhat arbitrary but are unlikely to lead to 
errors exceeding a factor of three. 

The derived rates of energy loss for electrons with’ energies 
ranging from 10 keV to 20 eV are given in Table 4 for air of unit 
density and the derived rates corresponding to the two model 
atmospheres are shown in Figures 8 and S. 

3.2 Excitation to Metastable States 

At energies of about 20 eV, energy loss through excitation of 
the metastable levels of atomic and molecular oxygen and of molecular 
nitrogen are comparable to the energy loss through allowed transitions, 
becoming of increasing importance as the energy decreases. Many 
processes contribute to the energy loss and very few of the relevant 
cross section data are available. A rough analysis suggests that the 
mean rate of energy loss in the atmosphere decreases slowly to a value 
of the order of 10 ^ n eV cm ^ at an energy of about 7 eV. Below 
this energy no significant loss is to be expected from excitation of 
electronic levels of and the rate of -energy loss at altitudes below 
300 km will decrease sharply as the energy falls below about 7 eV. 

It follows from the cross section calculations of Seaton (1956) 
that between 7 eV and 3 eV excitations to the metastable X D and 
states of atomic oxygen contribute about 5 x 10 ^ n (0) eV cm 1 to the 
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rate of energy loss. The contributions of dissociative attachment of 
molecular oxygen (cf. Thompson 1959) and of excitation to excited 
electronic levels of molecular oxygen (Schulz and Dowell 1962) are 
small in comparison and can be ignored at altitudes above 120 km. 

The derived rates of energy loss for electrons with ’energies 
between 20 eV and 5 eV for the two model atmospheres are shown in 
Figures 10 and 11. 

3.3 Vibrational Excitation of 

At energies below about 3.5 eV, energy loss through excitation 
of the vibrational levels of molecular nitrogen becomes very , 
important. The total cross sections have been measured by Haas (1957) 
and the relative cross sections for excitation to the individual 
vibrational levels with quantum number v between 2 and 8 by Schulz (1962). 
The corresponding rate of energy loss for unit number density of 
nitrogen is illustrated in Figure 12. It includes an arbitrary 
estimate of the contribution from excitation to the v = 1 level, the 
cross section for which has not been measured. 

The rates of energy loss appropriate to the two model atmospheres 
are given in Figures 13 and 14 which include also the contribution of 
excitation to the metastable levels of atomic oxygen. 

The fate of the energy stored as vibrational excitation is 
uncertain. It may be transferred to the ambient electron gas by 
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.5 eV 



Figure 13. Rates of energy loss to neutral particles (solid lines) and 
to electrons (dashed lines) for electrons with energies 
between 1.5 eV and 4 eV moving in the model atmosphere with 
T = 1000°K . 




Figure 14. Rates of energy loss to neutral particles and to electrons 
for electrons with energies between 1.5 eV and 4 eV moving 
in the model atmosphere with T (oo) = 2000°K. 



superelastic collisions 


e + N 2 (v) -+ e + (v 1 < v) (14) 

or it may be degraded into kinetic energy of the neutral * particles , a possible 
reaction being atom-atom interchange with atomic nitrogen (Bates 1955, 
Dalgarno 1963) 


N + (v) (v* < v) + N . (15) 

The rate coefficient of (14) may be of the order of 10 ^ cm^ sec 
but the rate of (15) is unknown except for transitions in which v’ = 0. 
We shall investigate the consequences of both possibilities. 

3.4 Rotational Excitation of 

Vibrational excitation is rather inefficient at energies below 

1.5 eV and the measurements of Schulz and Dowell (1962) suggest that 

1 + 1 

excitation of the 2 and A states of molecular oxygen are very 

§ § 

slow processes. As the energy approaches thermal energies, excitation 
of the rotational levels of is probably the most efficient loss 
mechanism, the small quadrupole moment of 0^ causing excitation of the 
rotational levels of 0^ to be comparatively slow. The rate of energy 
loss by rotational excitation of can be obtained from the calcula- 
tions of Dalgarno and Moffett (1962a) and it is about 
-20 2 -1 

5 x 10 n (N^) eV cm sec The corresponding rates of energy loss 
for the two model atmospheres are shown in Figures 15 and 16. 
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Figure 15. Rates of energy loss through rotational excitation of N 
(solid lines) and through elastic collisions (dashed 
lines) with the ambient electrons for T (oo) = 1000°K) . 




Figure 16. Ratos of energy loss through rotational excitation o 
(solid linos) and through clastic collision;; (dashed 
lines) -with the ambient electrons for T (oo) = ?000°K 



Z.5 Elastic Scattering by Neutral Particles and 'Positive Ions 

The rate of energy loss through elastic scattering by neutral 
particles is given by 

- ;r = tt n E Q (16) 

dx M m 

where m is the mass of the electron, M is the mass of the neutral 
particle and Q is the cross section for momentum transfer. We shall 
not present values of (16) since except at energies very close to 
thermal energies energy loss through rotational excitation is much 
more efficient. At high altitudes where the concentration of atomic 
oxygen exceeds that of molecular nitrogen, this is no longer the case 
but there, loss through elastic scattering by the positive ions is a 
more efficient mechanism. It is clear, however, from the mass factor 
in (16) that energy loss through elastic scattering by the positive ions 
is negligible compared to that through elastic scattering by the 
ambient electrons. 

Although elastic scattering by heavy particles can be ignored in a 
consideration of the slowing down of the fast photoelectrons, it 
provides an important mechanism for cooling the heated electron gas and 
we shall return to it in Section 4. 

3.6 Elastic Scattering by Ambient Electrons 

In discussing the relative efficiencies of various processes 
which slow down energetic electrons, we have ignored the elastic 
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scattering by the ambient electrons, an important process for the 
selective heating of the ambient electron gas (Hanson and Johnson 
1961). 

For the mean rate of energy loss through elastic collisions with 
the ambient electrons, we adopt a formula of Butler and Buckingham 
(1962) 


_1_ dE 

n dx 
e 


5.01 x 10 
Vie/ 1000 ' 


-ZZ.F ( 3 /41 E > ev cm 2 

VE V Te/ 10C0 


(17) 


where 


1 2 2 
F(x) = — • J ex P ^“ x ^ dx “ 2 exp ( -x ). 


( 18 ) 


o 

The value of Zn \ appropriate to the ionospheres of Table 1 is about 
15; it is unnecessary for our purposes to take account of its detailed 
variation with electron temperature and electron density. 

Values of (17) are given in Table 5 for a range of incident 
energies and for several electron temperatures. For energies greater 
than 3 eV, (17) may be replaced by the asymptotic formula 


i dE 1.95 x 10 

n dx E 


eV cm 


2 


e 

The energy loss rates through elastic collisions with the ambient 
electrons are compared in Figure S-16 with the energy loss rates 
corresponding to collisions with the neutral particles. 


( 19 ) 
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is taken to be 15. 


3.7 Critical Altitudes and Energies 


The results of the comparisons of the efficiencies of energy 
loss to the neutral particles and to the electrons are summarized in 
Table 6 which gives for various electron energies the critical 
altitudes for the two model atmospheres above which loss to the 
electrons is the more probable and below which loss to the neutral 
particles is the more probable, it being assumed that vibrational 
energy of is converted to thermal energy of the neutral particle 
gas. If we assume that vibrational energy of "S is converted to 
thermal energy of the electron gas, the critical altitudes are 
unchanged for E > 5 eV but do not exist for lower energies, all of 
the electron kinetic energy below 5 eV being transferred to the 
ambient electrons. 

We may summarize the results also by introducing critical 

energies E which are such that for E > E , E is transferred to the 
c c c 

electron gas and for E < E^, E is transferred to the electron gas. 

The critical energies are shown in Figures 17 and 18 for the two 
model atmospheres corresponding to the alternative assumptions that 
vibrational energy of is converted to thermal energy of the 
neutral particle gas and that it is converted to thermal energy of the 
electron gas . 
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able 6 


Critical Altitudes 


T (co) = 

1 , 000°K 

T (cc) = 

2 , 000°K 

Energy (eV) 

h (km) 
c 

Energy (eV) 

h (km) 
c 

40 

500 

18 

500 

20 

370 

14 

408 

10 

286 

10 

332 

8 

260 

8 

276 

7 

236 

7 

252 

6 

178 

6 

176 

5 

170 

5 

170 

4 

220 

4 

230 

3 

224 

3 

240 

2 

204 

2 

210 

1.5 

152 

1.5 

150 
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500 



Figure 18. Critical energies for the model atmosphere with T (V) =: 2000°K, 
assuming that vibrational energy of N 9 is converted to thermal 
energy of (a) the neutral particles and (b) the electron o a s. 




4. Heat Flux and Cooling Processes 


Combining the predicted initial energy distribution of the 
photoelectrons with the critical altitudes and energies yields the 
distribution of energy which is transferred to the ambient electrons. 
The results are shown in Figur.es 19 and 20, there being four 
distributions for each atmosphere depending upon the ultimate fate 
of the energies initially stored as metastable energy and as 
vibrational energy . 

We assume that elastic collisions between the ambient electrons 
are sufficiently frequent that the electron velocity distribution is 
Maxwellian, characterized by an electron temperature Te. The 
electron gas cools by collisions with the neutral particles and 
positive ions, which we assume have a common temperature Tg. For 
the energy loss in elastic collisions with atomic oxygen we adopt the 
expression 


(^-r~) = ~ 10' 14 n (0) Te^ (Te - Tg). 

dt o & 


( 20 ) 


which corresponds to a constant collision cross section of 6 x 10^ 

2 

cm (Rothe, Marino, Neynaber and Trujillo, 1962). For the energy loss 
in collisions with molecular nitrogen we use 


d Te 


/ . 


-16 


dt 'N 


>„ = - (7.6 x 10 Te + ( 


G v 


■) 


n(N^) rot 


n(N 2 )(Te - Tg) , (21) 
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the elastic loss following from the measurement of Pack and Phelps 

(1961) and the rotational loss parameter (Gv/n) from an extension 

rot 

of the calculations of Dalgarno and Moffett (1962b). For loss to 
molecular oxygen, we use an elastic contribution only 


( 77 )- = - 4.7 x 10 ' 16 n (0.) Te (Te - Tg) 

at 2. 

(cf. Dalgarno 1961). The energy loss in elastic collisions with 
positive ions is significant above an altitude of 200 km. At such 
altitudes we may suppose the positive ions are mainly 0 and the 
corresponding rate of energy loss is given by 


( 22 ) 


dT . _ % <Ie ~ Tg) 

dt + 268 T <? U 

(cf. Spitzer 1956, Hanson and Johnson 1961 ). If Q is the heat in 
-3 

ergs cm transferred to the electron gas in unit time, then in 
equi librium 


(23) 




3 kn 


Q = - 


(Hi) + (il®) + (l-l®) 

dt 0 dt N t 2 dt 0 2 


* <^> . i 


dt + 


(24) 


Because it ignores the contributions from vibrational excitation, 
the adopted cooling rate is probably too slow at high temperatures. 
However, the uncertainties in the cooling rates are smaller than the 
uncertainties in the predicted heat. fluxes. 
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5. Electron Temperatures 


It is instructive to consider first the equation which results 
when cooling occurs only by elastic collisions with positive ions. 

Q is measured in eV cm sec , the simplified equation for the electron 
temperature is 


Q = 5.51 x 10 



(Te - Tg) 



(25) 


As Hanson and Johnson (1961) have remarked, (25) has two solutions for 
Q less than a critical value Qci but for Q greater than Qci (25) has 
no solutions and Te is not limited by energy transfer to the positive 
ions. The critical heat input is given by 


Qci = 2 x 10 ^ n “ Tg 2 eV cm J sec x (26) 

6 

and the corresponding equilibrium electron temperature is given by 

Te = 3 Tg (27) 

The modifications caused by the inclusion -of cooling by energy 
transfers to the neutral particles depend upon the altitude. Below 
about 200 km, cooling in collisions with neutral particles is more 
efficient at all temperatures than cooling in collisions with positive 
ions and there is only one solution to (24). At greater altitudes, the 
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position is more complicated. Figure 21 shows the right-hand side or 

Equation (24) as a function of Te at an altitude of 400 km in the 

atmosphere with an exospheric temperature of 1000°K. It is apparent 

4 -3 -1 

that for heat inputs Q greater than 1.3 x 10 eV cm sec and less 
3 -3 -I 

than 7.8 x 10 eV cm sec there is only one solution of (24) but 

that for.O lying between these values there are three solutions. Gf 

the three solutions, the smallest and the largest -are stable in the 

sense that after a small disturbance the departure of Te from its 

equilibrium value will decay exponentially with increasing time and 

the remaining solution is unstable. 

It is also apparent from Figure 21 that as Q passes through the 
4 -3 - 1 

value 1.3 x 10 eV cm sec from below the equilibrium electron 

temperature increases discontinuously to a very large value. The 

3 

opposite behaviour occurs as Q passes through the value 7.8 x 10 eV 
3 -1 r 

cm sec from above. 

Solutions of Equation (24) are given in Table 7 for a wide range 
of heat inputs, together with the critical heat inputs Qci. The 
equilibrium electron temperature increases very slowly as Q increases 
through several orders of magnitude until Q exceeds Qci when an abrupt 
increase in Te occurs at high altitudes. The very large values of Te 

obtained at high altitudes for Q greater than Qci do not have any 

quantitative significance (since at such values the adopted cooling 
rate is much too small and heat conduction in the electron gas is an 

important effect) but they show that for Q in the neighborhood of Qci 
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Table 7 


Equilibrium Electron Temperatures for a Range of Heat Inputs Q 

T (oo) = 1000°K 



0 

Altitude 

(km) 

120 

150 

200 

300 

400 

450 

(eV 

-3 -1. 

cm sec ) 

Tg (°K) 

380 

763 

952 

998 

1,000 

1,000 


10 1 

Te (°K) 

380 

763 

953 

998 

1,000 

1,000 


io 2 


381 

767 

960 

1,000 

1,000 

1,000 


io 3 


385 

801 

1,040 

1,020 

1,030 

1,050 


io 4 


427 

1,160 

2,390 

1,380 

1,680 

132,000 


IO 5 


847 

4,190 

12,000 

43,700 

223,000 

— 

T (cc 

0 = 2000°K 

Tg (°K) 

380 

880 

1,390 

1,810 

1,940 

1,570 


io 1 

Te (°K) 

380 

880 

1,390 

1,810 

1,940 

1,970 


io 2 


381 

884 

1,400 

1,820 

1,950 

1,930 


io 3 


385 

924 

1,510 

1,870 

2,030 

2 , IGu 


io 4 


427 

1,330 

2,920 

2,800 

6,390 

22,600 


10 5 


847 

4,590 

11,900 

24,600 

64,100 

103,000 
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a large increase in Te may result from a small increase in Q. 

Corresponding to each of the predicted hear inputs, shown in 
Figures 19 and 20, there is only one possible equilibrium electron 
temperature and it is stable. The distributions of Te • with altitude 
are shown in Figures 22 and 23 from which it appears that the 
difference between the electron temperature and the heavy particle 
temperature at noon in a quiet ionosphere may attain a maximum value 
of between 800°K and 1500°K (depending upon the role of metastable 
and vibrational energy) at an altitude near 250 km and that it vanishes 
below about 120 km and above about 400 km. Very large values of Te 
do not occur because the predicted heat inputs exceed the critical 
values Qci only at altitudes below about 250 km where cooling by energy 
transfer to neutral particles is at least comparable in efficiency at 
all temperatures to cooling by energy transfer to positive ions. 

The predicted values are sufficiently large to suggest that the 
heat inputs at high altitudes might exceed the critical values near 
dawn and during severely disturbed conditions with the consequence 
that very large values of Te may occur then, although the 
heating effect of the increase in Q will be reduced by the associated 
increase in the ambient number density. Simultaneous measurements 
of Te and at dawn would be of great interest. 
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6. Discussion 


The observational evidence concerning electron temperatures in the 
ionosphere has been reviewed recently by Bauer and Bourdeau (1962). 

They conclude that except at sunrise the electron and neutral particle 
temperatures in a quiet ionosphere are equal up to altitudes of 200 km 
and at altitudes above the F2 peak. In the lower F region, the 
difference between the electron temperature and the heavy particle 
temperature attains a maximum value of about 1000°K at an altitude 
near 250 km. 

In view of the several uncertainties in the theoretical analysis, 
the agreement between the theoretical predictions and the observations 
is very satisfactory. Indeed the observed magnitude of the maximum 
temperature difference is closely equal to that predicted for the case 
when vibrational energy of molecular nitrogen is not transferred to 
the electron gas. The predicted location of the maximum is about 30 km 
lower than that observed but this small discrepancy can plausibly be 
attributed to an underestimate of the cooling rate, especially in view 
of recent work by Phelps (1962, private communication) indicating that 
the cross section for low energy excitation of the first vibrational 
level of molecular nitrogen is quite large. Thus the observational 
data on electron temperatures in a quiet ionosphere can be explained 
as due to the direct action of solar ultraviolet radiation and it is 
unnecessary to postulate the existence of electric fields- 
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At sunrise and during disturbed conditions the observational 
results indicate greater departures from temperature equilibrium 
(cf. Bauer and Bourdeau 1962) and work is in progress to analyze these 
phenomena. 


Geophysics Corporation of America 
Bedford, Massachusetts 
November 1963 
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